Ballistic thermal rectification 
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We study ballistic thermal transport in three-terminal atomic nano-junctions by the nonequilib- 
rium Green's function method. We find that there is ballistic thermal rectification in asymmetric 
three-terminal structures because of the incoherent phonon scattering from the control terminal. 
With spin-phonon interaction, we also find the ballistic thermal rectification even in symmetric 
three-terminal paramagnetic structures. 
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Introduction Phononics, the study of information pro- 
cessing and controlling of heat flow by phonons, is an 
emerging new field that is attracting increasing attention 
[ij. Specifically, researchers have recently modeled and 
built thermal rectifiers Q, thermal transistors Q, ther- 
mal logical gates and thermal memory which are 
the basic components of functional thermal devices. The 
most fundamental phononic component is the thermal 
rectifier - a device that allows larger conduction in one 
direction than in the opposite direction when it is driven 
far enough from equilibrium. The effect of thermal recti- 
fication is well known to be realized by combining the sys- 
tem inherent anharmonicity with structural asymmetry 
. Whether the rectification can happen in harmonic 
systems, that is, ballistic rectification, is still unknown, 
although one recent paper [l^] discussed it, in which the 
authors get different heat conductances by exchanging 
the heat baths with different structures connected to an 
asymmetric center part, thus the different conductances 
come from the totally different systems. As we know, the 
thermal transport in nanoscale materials, which is very 
promising for thermal devices, can be regarded as bal- 
listic because of their small sizes in comparison with the 
phonon mean free path. Therefore, it is highly desirable 
to investigate whether the ballistic thermal rectification 
can be realized in harmonic systems, and to explore the 
necessary conditions for thermal rectification. 

The ballistic thermal transport in two-terminal junc- 
tions can be described by the Landauer formula. Since 
the temperatures enter only through the Bose distribu- 
tion, it is obvious that if we reverse the heat bath tem- 
peratures, the heat flux only changes sign, and no rec- 
tification is expected. How about the ballistic thermal 
transport in multiple-terminal junctions? The theory 
for multiple-terminal electric transport was proposed as 
the Landauer-Biittiker conductance formul a |llMl3l| , and 
was applied to thermal transport recently IJ- la|. F rom 
the electronic transport in three-terminal system (l3| , we 
know that the third terminal can introduce incoherence 
or phase breaking to the transport. So it is our inter- 
est to investigate whether a multiple-terminal junction 
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FIG. 1: (Color online) The three-terminal junction setup to 
study the ballistic thermal transport. The left and right leads 
have temperatures Tl and Tr, the control terminal lead is 
adjusted to be Tc or so that the heat flux from this lead 
is zero in the forward (Tl = T+, Th = T_) or backward {Tl = 
r_, Tr = r+) process. T+ and T- are the temperatures of 
the hot and cold baths, respectively. 



is a proper option for ballistic thermal transport, that 
is, whether the incoherence through the third terminal 
can induce rectification effect. We will take the nonequi- 
librium Green's function (NEGF) approach d 
NEGF is widely applied to electronic and thermal trans- 
port, and is successful to study the spin Hall effect and 
phonon Hall effect injunctions [l6l[2g|. 

Model and Method We consider the ballistic thermal 
transport in a three-terminal nano-junction as shown in 
Fig. [U where a two-dimensional lattice sample, which 
is a honeycomb lattice, is connected with three ideal 
semi-infinite leads. The masses are coupled through 
nearest neighbors by elastic springs (with longitudinal 
and transverse force constants). We denote the center 
lattice as Nji x Nc, Nn, Nc correspond to the num- 
ber of rows and columns, respectively. The external 
magnetic field can be perpendicularly applied to this 
part. We use Ncl to denote the number of columns 
of the control lead and Nqd to denote the number of 
columns deviate from the middle of the center part; if 
NcD = 0, the whole setup is symmetric. In Fig. [l] 
Nr = 9, Nc = 8, Ncl = 4,, Ncd ^ -2. The Landauer- 
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Biittiker formula can be expressed as, 
duj 



J a 
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flUJTpa{'^)[n{Ta) - n{Tp)\. (1) 



Here, Tp^ is the transmission coefficient from the ath 
bath to the /3th bath; and niTa) = [e^^/^BT^ _ jg 
the Bose distribution with being the temperature of 
the ath heat bath. We set Ti ^ 1 and fc^ = 1 in the 
foUowing calculation. Therefore, in the forward process, 
Tl = T+ , Tff = r_ , we obtain 

j'^{TRL{oj)[n{T+) - n(r_)]+rcLH[n(T+)-n(rc)]} 

(2) 

Jc= J ^{tlcHHTc) - niT+)]+raciio)[niTc)-n{T-)]} 

Similarly we can obtain the heat fluxes J_ and in 
the backward process = T_, Tr =1+. From the 
equations of Jc = and J, 



c 



we can obtain the 
temperatures of the control bath Tc and T^; inserting 
them to the formulae of J+ and J_, by the definition of 
rectification as 



R — (J+ — J_)/max{J+, J_}, 



(4) 



we can calculate the rectification of this model. If the 
system is in the linear response regime or in the classic 
limit, the heat flux from the ath lead can be expressed 
as J„ = J2 ^pc,{T^ - Tp). we set r+ = Tq + A, T_ = 

To - A, Tc = To + (5, T'^Tq + 5\ then we obtain 



6=-5' = 



J+ = J_ = 2A aRL 



(TLC + CTRC 



(5) 



(6) 



So there is no rectification in the linear response regime or 
in the classic limit. In order to get thermal rectification, 
we should consider the quantum regime out of linear re- 
sponse, and the key work is to compute the transmission 
coefficients among the heat baths. 

In this paper, in addition to the structural asymmetry, 
we also can introduce the same spin-phonon interaction 
as in Refs. 1^ 22-24| in order to break a time-reversal 
symmetry. The Hamiltonian of our model can be written 
as 



H 



E 



E 



C/oMPo, 



(7) 



I3=L,R,C 



where = I {P^Pa + UlKJJa). The superscript T 
denotes matrix transpose. Here, 0, L, R, C correspond 
to the center region, left, right, and control leads, re- 
spectively. Ua is a column vector for mass reduced dis- 
placements in region a. Pa is the associated conjugate 
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FIG. 2: (Color online) Rectification as a function of relative 
temperature difference of the hot and cold heat baths. The 
parameter of the setup is Nr — 9, Nc ~ 16, Ncl = 2. The 
temperature of the heat bath are T+ = ro(l -I- A) and T- = 
ro(l — A), where To = 0.2 is the mean temperature. The solid 
square, solid circle, solid triangle, diamond, hollow triangle, 
hollow circle, hollow square correspond to Ncd = —7, -5, -3, 
0, 3, 5 and 7, respectively. 



momentum vector, and is the force constant matrix, 
^/3,o = (^,^)"^ is the coupling matrix between the /3th 
lead and the central region. A is a block diagonal matrix 

with diagonal elements ^ ^ ^ . /i is a model parame- 
ter which is supposed to be proportional to the magnetic 
field. In the first part of this paper, we set /i = 0; it is a 
standard harmonic phononic system. 

The retarded Green's function for the central region in 
frequency domain is 16 1 



G"'M= {uj + irif - Ko-T.''[uj]- + 2iojA 



(8) 



Here, S'' = J2 ^/3' ^^'^ — ^,^5^^/3,0 is the self- 

I3=L,C,R 

energy due to interaction with the heat baths, gj^ = [{uj + 
ir])'^ — Ki3]~^. T] is an infinitesimal real positive quantity. 
The surface Green's functions of the leads g"^ are obtained 
following the algorithms of Ref. [isj. The transmission 
coefficient is 

T^„M = Tr(G'T^G'T„), (9) 
where, F„ = ^(E^M - E^M), G« = (G'^)t, and - 

Results and Discussions Firstly, we consider the ballis- 
tic thermal transport in an asymmetric structure without 
an external magnetic field. In the following simulation, 
we set the longitudinal spring constant — 1.0, and 
the transverse one kr = 0.25. If the control lead is con- 
nected to the middle of upper edge of the center, that 
is, Ncd = 0, the forward process and backward one are 
exactly the same; no rectification will be expected, as 
shown in Fig. [2] (the diamond symbols). We will obtain 
the same result of no rectification if we replace the cen- 
ter honeycomb lattice with the square lattice same as 
the leads, when the whole system is symmetric wherever 
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FIG. 3: (Color online) The diflterence of transmission coeffi- 
cients: tlc — trc, as a function of frequency. The parameter 
of the setup is Nr = 9, iVc = 16, Ncl = 2. The solid, dot 
curves correspond to Ncd ~ —3 and Ncd ~ —7, respectively. 



we put the control lead. If the control lead moves avifay 
from the center, the rectification effect appears. When 
the lead is moved the same distance to the left or right, 
the rectification coefficient has the same magnitude but 
opposite sign, which is because that the two cases only 
exchange the value of J+ and J_. If the distance be- 
tween the control lead and the middle of the center part 
is longer, the rectification effect is larger. In Fig. [5J we 
can see that the case of Ncd = ±7, when the control lead 
is next to left or right lead, has biggest rectification. The 
rectification increases with the temperature difference at 
far- from- linear-response regime. 

From the above formulae, we find if the transmissions 
Tlc and trc are linearly dependent (here Tap = Tpa be- 
cause of time reversal symmetry), the rectification will 
be zero. Figure [3] shows that the transmission coeffi- 
cients Tlc and trc are not linearly dependent, so that 
the rectification is nonzero. In some frequency domains, 
Tlc > TRC] and in other frequency domains, tlc < tuc- 
With the distance increasing, that is, from Ncd = —3 
to Ncd = — 7, the difference of tlc and Tnc enlarges 
in most part of the whole frequency domain; so that the 
rectification coefficient increases. Because of scattering 
from the third bath - the control lead, the total thermal 
transport from the left lead to the right one is partially 
incoherent. This phonon incoherence induces rectifica- 
tion. In the whole system the Hamiltonian is quadratic, 
that is, there is no nonlinearity or anharmonicity, but we 
still can obtain rectification. Therefore, the phonon in- 
coherence, which can be induced by either nonlinearity 
or scattering lead, is the necessary condition for thermal 
rectification. 

The control lead acts as a scattering source, which 
makes the phonon transport incoherent, so that the rec- 
tification comes out. However, the width of the con- 
trol lead does not quite affect the whole thermal trans- 
port, which is shown in Fig. IH^a). We make the control 
lead next to the left lead, and find that the rectification 
changes little when we increase the width of the lead. If 
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FIG. 4: (Color online) (a) Thermal rectification as func- 
tion of relative temperature difference A for different width 
of control lead, at To = 0.2; (b) Thermal rectification as func- 
tion of mean temperature for different relative temperature 
difference. Ncl = 2, Ncd ~ —7. For both (a) and (b), 
Nr = 9,Nc = 16, T+ = 70(1 + A) and r_ = ro(l - A). 



the width of control lead increases further, the rectifi- 
cation decreases because the asymmetry decreases. Fig- 
ure Hljb) shows the rectification dependence on temper- 
ature, and reproduce the reversal of rectification found 
in Ref. I9|. At a low temperature, the contribution to 
thermal transport only comes from the low frequency 
phonons; if the temperature increases, more high fre- 
quency phonons will contribute to the heat transport. 
From Fig. |3l the relations between transmissions tlc and 
TLic in low frequency domain and high frequency domain 
are opposite, so that the rectification reverses with the 
temperature increasing. When the temperature increases 
further, the system will go to the classic limit, the recti- 
fication disappears. 

From the previous work fi], 0, [HI on thermal rectifica- 
tion, we know that in order to get rectification, we need 
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FIG. 5: (Color online) (a) Thermal rectification as function of 
relative temperature difference A for different external mag- 
netic fields, (b) Thermal rectification as a function of mag- 
netic field h. For both (a) and (b): Nr = 9, Nc = 16, Ncl = 
2, Ncd = 0. r+ = ro(l + A), T- = ro(l - A). ro=0.2. 
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the structural asymmetry. However, in the nanoscale 
rectifier, it is not easy to control the structural asym- 
metry or not easily distinguish the rectification direction 
by the structural asymmetry. Is there any other means 
to introduce asymmetry to induce rectification? From 
the study of phonon Hall effect 1^ 22-24|, it is known 
that the magnetic field can infiuence the thermal trans- 
port by the spin-phonon interaction. Thus the magnetic 
field can break the symmetry of the phonon transport. 
We apply an external magnetic field perpendicular to the 
center part of a symmetric structure to study the ballistic 
thermal transport, the results are shown in Fig. [5l The 
thermal rectification effect as a function of the tempera- 
ture difference is shown in Fig.[5ja). R increases with the 
temperature difference, and can be about 3% if A = 0.8 
and h = 0.3 at Tq = 0.2. Figure[5fb) shows that the rec- 
tification can monotonically increase with the external 
magnetic field in the range of /i = ^ 0.3. The ap- 
plied magnetic field breaks the symmetry of the phonon 
transport system through the spin-phonon interaction, so 
the transmission coefficient from the control lead to left 
one is different from that to right one. Figure IHl shows 
that the two transmission coefficients are not linearly de- 
pendent; and with the increasing of magnetic field, the 
difference of these two transmissions will enlarge, which 
induce bigger rectification effect. Here all quantities are 
dimensionless, if we use the parameters for real materials, 
the rectification coefficient may be different, but would 
be likely in the range of few per cents. 

Although all our values for rectification are small, if 
we can devise a chain of such ballistic rectifiers, the to- 
tal rectification can be as large as we need. Because 
for most nanoscale materials the thermal transport is 
ballistic and temperature can be applied far from lin- 
ear response regime, our prediction can be tested by ex- 
periments. The asymmetric quasi-two-dimensional nano- 
structure of any material can have ballistic thermal rec- 
tification, such as a graphene sheet. Or the symmetric 
structure with spin-obit interaction, such as paramag- 
netic dielectrics, can also have ballistic rectification ap- 
plied an perpendicular magnetic field. 

Conclusions Using the nonequilibrium Green's func- 
tion method, we have studied ballistic thermal transport 
in three-terminal atomic nano-junctions. Adjusting the 
temperature of the control lead in order to make the heat 
flux from this lead be zero, we can calculate the thermal 
rectification effect. In the quantum regime out of linear 
response, there is ballistic thermal rectification in asym- 
metric three-terminal structures because of incoherent 
phonon scattering from the control terminal. Through 
spin-phonon interaction, we also find the ballistic thermal 
rectification in symmetric three-terminal paramagnetic 
structures applied an external magnetic field. Therefore, 
not the nonlinearity, but the phonon incoherence, which 
can be induced by nonlinearity or scattering boundary 
or scattering lead, is the necessary condition for thermal 
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FIG. 6: (Color online) The difference of transmission co- 
efficients: tlc — trc, as a function of frequency for differ- 
ent applied magnetic fields. The parameter of the setup is 
Nr = 9, Nc ^ 16, NcL = 2, NcD = 0. The solid, dot curves 
correspond to h = 0.1 and h = 0.3, respectively. 



rectification. Another necessary condition is asymmetry, 
not necessarily being structural asymmetry, which can be 
introduced by an applied external magnetic field through 
the spin-phonon interaction. 
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